PROCEEDINGS 

™ SCIENCE 



Long distance contribution to Kl-Ks mass difference 



o 

(N 

> 
O 



i3 



> 

in 

On 



Jianglei Yu"^ 

Department of Physics, Columbia University, New York, NY 10025, USA 
E-mail: jy2 37 9@columbia.edu 

We present a method to non-perturbatively determine the long-distance contribution to the Ki- 
Ks mass difference. The calculation is performed on 2+1 flavor, domain wall fermion, 16^ x 32 
configurations with a 421 MeV pion mass and a kaon mass of 559 MeV . We include only current- 
current operators and drop all disconnected diagrams in the calculation. The largest contribution 
comes from quadratically divergent, short distance lattice artifacts. This quadratic divergence is 
eliminated through the GIM mechanism by introducing a valence charm quark. The remaining 
short distance effects are then removed by using RJ/MOM technique which allows their exact 
replacement by the physical short distance part. 



X 
S3 



XXIX International Symposium on Lattice Field Theory 
July 10-16 2011 

Squaw Valley, Lake Tahoe, California 



* Speaker. 



© Copyright owned by the author{s) under the terms of the Creative Coinmons Attribution-NonCoinmercial-ShureAHke Licence. 



http: //pos. sissa. it/ 



Long distance contribution to K^-Ks mass difference 



Jianglei Yu 



1. Introduction 

The extremely small K^-Ks mass difference has been measured accurately decades ago. The 
origin of this difference is the K^-K^^ mixing via second order weak interactions. Convention- 
ally, the mass difference is separated into short-distance part and long-distance part. While the 
short-distance contribution has been calculated to next-to-leading order [1], the long-distance con- 
tribution could only be determined non-perturbatively, which contributes around 30% to the mass 
difference [2]. Norman Christ suggested a Lattice QCD method to compute the long-distance con- 
tribution [3]. This proceeding is the first numerical experiment about the new method. 



2. Second order correlator 

To compute Kl-Ks mass difference on a Euclidean space lattice, we could evaluate the time- 
integrated second-order product over a time interval , tt] : 

^=\tt< K\tf)Hw{t2)Hw{h)K''Hh) > (2.1) 

tl=ta t2=ta 

Here the Kaon is created at t, , the second-order weak Hamiltonian acts between time interval ,th\, 
and the outcome anti-Kaon is annihilated at tf. The ampUtude is represented schematically in 
Figure 1. 




Figure 1: Demonstration of the second order correlator in Equation 2.1. Here the two four quark opera- 
tors are integrated over the shadowed region. 

If we assuming that ta — ti and tf — ti, are large enough for the interpolating operators to project 
onto Kaon states, after inserting a complete set of intermediate states, Eq 2.1 becomes: 



|2^|2^-M^(r,-r,) 



^ <K \Hw\n><n\Hw\K^> ^^(m,-m„)t , 



1 



(2.2) 



{MK-E„r 

Here T = th — ta + l and Zk is the normalization factor of Kaon interpolating operator. We assume 
that there is no intermediate state degenerate with kaon in this expression, which is true in this 
work. The term proportional to T in Equation 2.2 gives the finite volume approximation to Kl-Ks 
mass difference. 



y <K"\Hw\n ><n\Hw\K^ > 

r MK-En 

The other terms in Equation 2.2 on T could be classified into three categories: i). The exponential 
decreasing term, if > Mk, these terms could be neglected when T is sufficiently large; ii). Ex- 
ponential increasing term, if E„ < Mk, these terms must be identified independently and subtracted 
from the result; iii) The term independent of T, which is trivial. 
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The full AS = 1 effective Hamiltonian consist of 7 independent four-quark operators [4], we 
include only the current-current operator Qi in this work. 

Qi = isada)v-AiupUj})v-A (2.4) 

The four different types of contractions are listed in Figure 2. In this work, we only compute type 
1 and type 2 contractions and drop type 3 and type 4 contraction. Type 3 contraction is dropped 
because it is disconnected diagram. We need to compute extra random source propagators for type 
4 contraction. So we also drop type 4 contraction in this first numerical experiment. 




Types 



Figure 2: Four types of contractions in the 4-point correlator, only type 1 and type 2 are included in the 
calculation, type 3 and type 4 are dropped. 

As we mentioned before, we must identified the exponential increasing terms in Eq 2.2. These 
terms comes from the intermediate states which are lighter than Kaon, such states are tt" state and 
vacuum state in this calculation. Since the disconnected diagrams are neglected, there will be no 
vacuum intermediate state. Then we must calculate < 7r°|2i|iC'° > and subtract the exponential 
increasing contribution from Equation 2.2. The contractions in this calculation are given in Figure 
3. We drop type 2 to be consistent with 4-point correlator calculation. 




Type 1 Type 2 

Figure 3: Two types of contractions in < 7i'^\Qi > 3-point correlator, type 2 is dropped to be consistent 
with 4-point correlator in Figure 2. 



3. Short distance effect 

The computation was performed on A'y = 2 + 1 flavor 16^ x 32 x 16 lattice with DWF, Iwasaki 
gauge action, a^^ = 1.73(3) Gev, 421 Mev pion mass and 559 Mev kaon mass. The two wall- 
source kaons ai^e located at time slice = and tf = 27. The two A5 = 1 operator act between 
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time slice [4,23]. We calculated 600 configurations separated by 10 Monte Carlo time units. The 
result is given in the first plot in Figure 4. The plot shows the integrated second order coiTclator 
as a function of integration time interval T. For each given integration time interval T, we calculate 
the integrate coiTclator in Equation 2.2 between all possible time interval [ta,ta + T — I] and take 
the average as final result. The two curves in the plot are the results before and after the subtraction 
of 7r° exponential term. The results have both long-distance part and short-distance part. The short 
distance part means that the two = 1 operator are close to each other on the lattice. We expect 
the short distance part to be quadratically divergent because of the up quark loop in Figure 1 . To 
get a detailed understanding of the short distance effect, we could introduce an artificial cut off, 
i.e., require the separation between two operators \x2 —xi\ > r. This cutoff will reduce the short 
distance effect and the long distance part will remain untouched. The other plot in Figure 4 show 
the results with cutoff radius 5. We can see that the amplitude of the result is reduced substantially 
after introcuding the cutoff. And the contribution from tt" intermediated state becomes visible. We 
could measure the mass difference at different cutoffs. From Equation 2.2, the mass difference is 
given by the coefficient of the linear term up to some factor when T is large. We choose to fit the 
slope of integrated coiTclator plot in the range T G [1 1 , 20] . The mass differences ai^e listed in Table 
1. We could do a naive inverse quadratic fit for the mass differences at different cutoffs. The result 
is in Table 1, which suggest the short distance effect is quadratically divergent. 



Table 1: Mass differences at different cutoff radius 



Cutoff Radius 1 2 


3 


4 


5 


M^K 0.3342(80) 0.1533(30) 


0.0796(17) 


0.0560(15) 


0.0455(14) 



4. Charm quark and GIM 

In order to remove the short distance, we introduce valence charm quark to the calculation. 
Then GIM mechanism will reduce the quadratic divergency into logarithmically. To implement 
this in the lattice calculation, we could replace all the up quark propagators in Figure 2 with the 
difference between up quark propagator and charm quark propagator. We use 5 difference valence 
charm quark mass ranged from 200Mev to lOOOMev. The integrated coiTclators after GIM subtrac- 
tion are plotted in Figure 6. The mass difference for difference valence charm quarks are listed in 
Table 2. 

5. Remain Short distance effect 

Even with the valence charm quark, there will still be some short distance lattice artifacts 
remain. The short distance part in Equation 2.1 could be given by: 

^SD = \t <^itf)C{^)mK'Hu) > 

^ t=ta 

= l-\ZK\^e-'^''^'f-'-'>C{^) < t^l^lK^ > T (5.1) 
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Figure 4: Integrated second order correlator as a function of integration time interval. Red and blue curve 
show the results before and after the subtraction of exponential term. First plot is the result without any 
cutoff. Second plot shows the result with cutoffs radius 5. 



Here T = th — ta + l, ^ = {sd)v-A{sd)v-A, is the conversion factor at a certain momentum 
scale jJ.. To identify C(/x) by using RI/SMOM technic, we could evaluate off-shell, four quark, 
amputated green function for the two diagrams in Figure 7 at some large external momentum 
scale jJ.. The external momentum satisfy = P2 = {pi — pi)^ = IJ-^ ■ Suppose the result for two 
diagrams are Gij^i and /^yjt/ respectively. Then we project the green functions into the desired 
gamma structure Pijti = {{^ — 7^)7n)ji{{^ — 7^)7^)ik- The conversion factor is given by Equation 
5.2. 

G{^) = GijkiPijki F{^) = FijkiPijti C(Ai) = ^ (5.2) 

In Figure 5, the left plot shows the dependence of C{lJ,) on the momentum scale jj. if we use 
1 Gev valence charm quark. As we expected, C(/x) will decrease while increasing n, because the 
difference between charm quark and up quark will decrease while the momentum scale get larger. 
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Figure 5: The mass differences at different cutoff radius, the blue cure is a naive two parameter fit. 



lOOOMev Charm Quark 



o 



u 


« B 






0.5 
-1 


# 1 
< 


1 

H 

n 

^ " " ™ □ . 


B D B S ^ ^ 


1.5 






a Original 








After subtraction 


-2 






2.5 








-3 






f 


3.5 






f 











5 10 15 20 25 

Integration Time Interval 

Figure 6: The integrated correlator after the inclusion of lOOOMev charm quark. Red and blue curve show 
the results before and after the subtraction of n'^ exponential term 




Figure 7: Off-shell, amputated four quark green functions, left diagram is for two AS' = 1 operators, right 
diagram is for one AS = 2 operator. 



6 



Long distance contribution to K^-Ks mass difference 



Jianglei Yu 





° 1Gev Charm quark 


I 


n - 

m 

m 


- I 





1.2 1.4 1.6 1.E 



n fi=2Gev 



Charm quark mass 



Figure 8: The left plot shows the conversion factor C(/i) as a function of momentum scale qith 1 Gev charm 
quark. The right plot shows C{iJ.) at different charm quark masses at fixed n = 2Gev. 



In the right plot, we fix momentum scale to be 2 Gev and plot C as a function of charm quark mass. 
When charm quark mass get smaller, the remain short distance effect becomes smaller. In Table 2, 
we show the remain short distance effect at different charm quark masses. We conclude that the 
remain short distance effects are so small that we could neglect them in this work. 



Table 2: Mass differences at different valence charm quark masses 



Mf(Mev) 


200 


400 


600 


800 


1000 


AMk 


0.0440(10) 


0.0455(12) 


0.0496(13) 


0.0556(14) 


0.0628(15) 


{AMk)sd 


6.2e-5 


2.4e-5 


6.2e-4 


0.0013 


0.0023 



6. Conclusion 

We perform a numerical study of the long distance part of Ki-Ks mass difference. The short 
distance part could be reduce by the inclusion of valence charm quark. The exponential increasing 
term could be identified and subtracted. The remain short distance effect could be computed by 
using RI/SMOM technic and removed from the results. 

The author thanks all his RBC/UKQCD collaborators for discussions and suggestions. Espe- 
cially thanks to Prof. Norman Chris for detailed instructions and discussion. 
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